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Abstract 

The first part of this paper develops a theory for the free energy of lyotropic 
polymer nematic liquid crystals. We use a continuum model with macroscopic 
elastic moduli for a polymer nematic phase. By evaluating the partition func- 
tion, considering only harmonic fluctuations, we derive an expression for the 
free energy of the system. We find that the configurational entropic part 
of the free energy enhances the effective repulsive interactions between the 
chains. This configurational contribution goes as the fourth root of the di- 
rect interactions. Enhancement originates from the coupling between bending 
fluctuations and the compressibility of the nematic array normal to the aver- 
age director. In the second part of the paper we use osmotic stress to measure 
the equation of state for DNA liquid crystals in 0.1M to 1M NaCl solutions. 
These measurements cover 5 orders of magnitude in DNA osmotic pressure. 
At high osmotic pressures the equation of state, dominated by exponentially 
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decaying hydration repulsion, is independent of the ionic strength. At lower 
pressures the equation of state is dominated by fluctuation enhanced elec- 
trostatic double layer repulsion. The measured equation of state for DNA 
fits well with our theory for all salt concentrations. We are able to extract 
the strength of the direct electrostatic double layer repulsion. This is a new 
and alternative way of measuring effective charge densities along semiflexible 
polyelectrolytes. 

PACS numbers: 87.15.Da, 61.30.-v, 64.30.+t 
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I. INTRODUCTION 



In biology, there exists a class of bulk materials that provides the cell's structural sta- 
bility and ensures its integrity in multi-cell environments. These materials are often made 
from biopolymers with varying intrinsic stiffness and helicity such as: actin, fibrin, col- 
lagen, polysaccharides, elastin, tubulin. In most cases biology controls elastic properties 
by controlled polymerization and depolymerization of the monomers and by enzymatically 
controlled crosslinking of polymer strands. Often biopolymers, because of their intrinsic 
stiffness, form liquid crystals at in vivo concentrations. Because these polymers are mechan- 
ically uniform, their phase behavior is predominantly determined by the volume fraction of 
the polymer in its solvent, rather than by temperature. Such materials are called lyotropic 
systems. In contrast, polymeric liquid crystalline materials used in industrial applications 
(like main-chain and side-chain polymer liquid crystals) are most often thermotropic, be- 
cause of their flexible backbones. Thermotropic means that temperature determines their 
phase behavior. 

Biopolymer materials are interesting for several reasons. First of all, there are many 
biomaterials made from or with biopolymers exhibiting mechanical properties unreached by 
most conventional synthetic materials, may the focus be strength, flexibility or a combination 
of both. Needless to say that little is understood about their design, structure and how they 
work. There is also the increasingly appreciated potential of biopolymers for basic material 
and condensed matter research fl[]. In most synthetic polymer systems it is extremely difficult 
to control properties of individual polymers, like their degree of polymerization, crosslinks 
and chemical uniformity. By using modern molecular biological and biochemical techniques 
it is possible to tailor biopolymers almost at will, using nature's own efficient polymerization 
machinery. So it is, for example, possible to prepare monodisperse DNA fragment solutions 
with lengths varying from a few nm to several fim using recombinant DNA methodology. 
This is not yet possible with any other polymer. 

DNA liquid crystals are also of immediate interest, because they model packing of DNA 
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in confined spaces, as in cell nuclei and viral phage heads (see e.g. [@-fi)- Moreover, it is 
believed that packing of DNA in chromatin plays an important role in gene- regulation ||. 
By understanding how much energy is needed to compact DNA we will gain insight into 
these processes as well. 

One of the questions in condensed matter physics is how to relate microscopic interactions 
to condensed state bulk properties, like elastic and dielectric constants, heat capacity and 
packing symmetries. In an experiment we can go the other way. The question is then: Can 
we infer microscopic interactions from macroscopic behavior? 

In this paper we explore the equation of state for DNA liquid crystals || using the 
osmotic stress method 0. It allows us to control all intensive variables of the system, 
namely the chemical potentials of both water and salt. According to the Gibbs-phase rule, 
this leaves us, after equilibration against solutions of known activities, with a single phase. 
The absense of phase coexistence is crucial for determining the structure of the different 
phases using x-ray scattering. At the same time the osmotic stress method provides us with 
information about the free energy of the system, which can be obtained by integrating the 
equation of state (osmotic pressure II versus DNA density). 

The ability to measure an equation of state and to map the phase diagram for free en- 
ergies made us realize that there was no polymer liquid crystalline theory for free energies. 
This is not surprising considering that typical liquid crystalline theories are continuum elas- 
tic theories expanding deformations around the symmetries of the system. Because these 
approaches are only good for long wavelength physics, it is very hard to use them to predict 
free energies. 

Despite these difficulties, we will show in the first part of this paper that in our partic- 
ular case with exponential interactions between the polymer chains, the free energy can be 
derived starting from a macroscopic theory. The microscopic interaction were incorporated 
into the continuum elastic moduli: the three Frank constants K\, K2, K% and the bulk com- 
pressibility perpendicular to the chains B. By evaluating the partition function, considering 
only harmonic fluctuations, we derive an expression for the free energy of the system. To 
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make the free energy finite we had to introduce a short wavelength cutoff. We argue that 
the form of the free energy is robust, meaning that the free energy behavior with respect to 
density does not strongly depend on the particular choice of the cutoff. 

We find that the fluctuation part of the free energy, going as the fourth root of the 
direct interactions, enhances the repulsive interactions between the chains. Enhancement 
originates from coupling between bending fluctuations and the compressibility of the nematic 
array normal to the average director. The results should be applicable to lyotropic polymer 
liquid crystals with fluid-like positional order like nematics or hexatics. We think that it even 
can be applied to chiral phases, like cholesterics and blue phases. In typical chiral phases 
the twist extends over several hundred intermolecular spacings, so that in these phases local 
parallel packing can still be used. 

We compare theoretical predictions with measurements of equations of state for DNA 
liquid crystals under different ionic condition over almost five orders of magnitude in os- 
motic pressure 0. By extracting the form of the fluctuation part of the measured free 
energy and assuming exponentially decaying repulsive interactions (screened electrostatic 
and hydration) between DNA molecules we have been able to extract the strength of the 
direct electrostatic double layer repulsion. This extraction creates a new and alternative 
way of directly measuring effective charge densities along semiflexibe polyelectrolytes. 



II. POLYMER NEMATIC THEORY 

Let us first consider the elastic free energy of a nematic: a three dimensional liquid 
with long-range orientational order and an average director n along the z-axis. Such phases 
are typically formed by solutions of rod-like or disc-like objects. There are three kinds of 
deformations in quadratic order of n with symmetry C^: splay, twist and bending. The 
corresponding elastic constants for these deformations are the Frank constants K\ (splay), 
K 2 (twist) and K 3 (bend). 

F N = U d 2 r ± dr z \K X (V • n) 2 + K 2 (n ■ (V x n)) 2 + K 3 (n x (V x n)) 2 j (1) 



For small deviations of the director field n(r) around its average orientation along the z 
axis n(r) ~ (Sn x (r), 5n y (r), 1), the free energy assumes the form 



T N = \ J d 2 r ± dz \Kt (V± • Sn) 2 + K 2 (V± x Sn) 2 + K 3 (djn) 



(2) 



For polymer nematics we now have to consider that the director field n(r) and the density 
of polymers in the (x,y)-plane p = po + Sp are coupled If the polymers were infinitely 

long and stiff the coupling is given by the continuity equation: 



d z 5p + po V± • Sn = 



(3) 



This constraint, however, is softened if the polymer has a finite length I or a finite 
persistence length £ p . On length scales larger than I or C p the polymer can either fill the 
voids with its own ends or fold back on itself [HI]. On these length scales the polymer 
nematic can splay without density change as illustrated in Fig.l. Following [|TT1 , p~2|1 this can 
be expressed by introducing G, a measure how effectively the constraint is enforced. Density 
changes are expanded to second order in density deviations 5p(r±,z) = p(r±,z) — p . B is 
the bulk modulus for compressions and dilations normal to the chains. The total free energy 
can be written 



T = MPo) + \ I d 2 ?±dz \b (|) 2 + G (d z 5p + p V ± ■ Sn) 



N 



(4) 



where G is given by 



G 
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with I never to exceed C p ||13|| . In the limit of finite polymer length G is also finite and can 
be obtained from the observation that d z 5p + po^7± ■ Sn equals the local difference between 
the number of chain heads and tails [[J . From here one derives that G is the concentration 
susceptibility for an ideal mixture of heads and tails thus 

k B T 



G 



(pH + Pt) ' 
6 



(6) 



where pn and pr are the average concentrations of heads and tails, with pjj, Pt = Pchain- The 
chain density on the other hand equals p c hain = ^f, wherefrom G = gp. The corresponding 
structure factor can be written as [|T^ 

5(g„g 2 ) = <\5p( q± , qz )f> = k B T - P ° 9 " * * (7) 

£gi + /c B T + q 2 z ) /C(q) 

where we defined 

K(q) = + (8) 
For long-fragment DNA the limit £ — > oo is appropriate, leading to the structure factor 



proposed by Selinger and Bruinsma [TJ 



S(q ± , q z ) = k B T 2 2 f q } 2 . (9) 

Xigl^ + K 3 q 4 z + Bq{ 

In order to calculate the contribution to the free energy due to fluctuations in nematic 
order we have to sum over all the density modes. From here on we will consider only 
fluctuations that are coupled to density changes. Because we are interested in how the free 
energy changes with density only those fluctuations will contribute. Twist deformations do 
not couple to density variations and can therefore be neclected. 

T = \ k B TJf d ^^ In [K x q\ql + K z q\ + Bq\) . (10) 

To evaluate the integral we take the partial derivative with respect to the compressibility B 

f)T 2 

_ _ I Ur.TV IT q± d q± d <iz q ± cni 

QB~ 2 B * ( 2 ^) 2 K^ql+K^+Bqi ■ y 11 ) 



The q z integral can be done straightforwardly and we remain with 



0Jr \k B Tj^l J — ^ — . (12) 



dB 2 " v ^i^/x l9 i+2 v ^i 

This integral depends essentially on the upper cutoff for q^ = q± max and we obtain 



where the function F(x) has been defined as 

r x u 3/2 du , , x !^ 5/2 )x < 1 

F(x) = / - = \ (^x~VT+x~(2x - 3) + 3 sinh- 1 ^) = { 5 . (14) 
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From here we obtain the two limiting forms of the free energy as 



k B T V J B 5/2 ^ BK 3 



f ~5x 23/% v k~ 3 q±max + ■ ■ ■ ;^™*« 2 y^2- ( 15 ) 

i6tt V 9±max + " ' '^ m ^ >2 y"^r- ( 16 ) 

Obviously the long-wavelength dependent physics is very complicated and depends crucially 
on the values of typical polymer length and the ratios of elastic constants. However it 
depends also on the q± cutoff. We have either to eliminate the cutoff by including higher 
order terms in the original Hamiltonian or to choose a meaningful cutoff. Higher order terms 
will capture the short-wavelength physics and remove the divergence. 

In the following paragraph we want to show how this can be done for the low density 



limit eq. [15|. We are aware that in a more thorough treatment one has to include all possible 
higher order terms, but for right now we only include terms that will make integral eq. |i~2| 
convergent. Letting B = B(q z , q±), K 3 = K 3 (q z , q±), K\ = K\(q z , q±) and taking into account 
the symmetry of the elastic nematic free energy, one can easily show that in order to make 
the integral in eq.|l^ converge in the low density limit eq. |15|, it is only necessary to expand 
B as 

= So (i + C 2 ?i + Al) ■ (17) 

All the other expansion terms either simply lead to a renormalization of K\ (presumed to 
be small anyhow) or they do not add up to the convergence of the free energy. 

This form immediately leads to a density- density correlation function in the (x,y)-plane 
which oscillates and decays exponentially |lBfl . Specifically the correlation function S(p,z) 



averaged over the length of the polymers assumes the form 



■^,) ■ (CV.J g^ l. (18) 
L Jo AuBqC 1 \r] { uf-ui I 



where the two complex zeros of 1 + u 2 + (2) u 4 = 0, and K (x) is the modified 

Bessel-function. The form eq.[18| of the transverse correlation function is a reasonable ansatz 
for a liquid structure factor. Evaluating the integral (eq. |TT| ) leads to the following expression 

k B T V Jb^ 5/2 r n i/2 d 



f = fo + -^ r \1 1 ^C b/2 — -ru (19) 



16v^tt V ^3 Jo 2 , f v \ 4 aV /4 



l+M 2 +(2)' 



u 4 



We essentially derived the same result as in eq. [L5|. The only difference is that q± ma x has 
been replaced by £ represents the correlation length in a liquid, which for a typical 
liquid goes as the distance between next neighors. Eq. [T^ suggests that instead of considering 
higher order terms we can simply choose a cutoff proportional to the Brillouin zone radius 
q±max — 5j where d is the effective separation between the polymers in the nematic phase. 
This is a physically meaningful and appropriate cutoff because the underlying macroscopic 
elastic model has, by definition, to break down at wavelengths comparable to the distance 
between molecules. 

Let us now see how our free energy scales for hard core repulsion between chains. First 



we take the limit of low density q± max <C 2^^r. For an elastic polymer K 3 has two additive 
terms, an intrinsic one and an interaction contribution. The first one stems from the elastic 
nature of the polymers themselves and has the form K3 ~ ksTCpPo, while the second term 
is the generic form valid for any steric interaction with a hard core a i.e. K 3 ~ k B T /a 
|T7| |. We now use the Helfrich self-consistent argument to evaluate the free energy due to 
elastic fluctuations in a hard core potential. It amounts to taking B = V^f?, where V is the 
volume of the system which can be taken as V = L x d 2 , L being the length of the sample. 
The crucial step is now to choose the cutoff which we set as q± max — 5- Inserting this into 
eq.pj^ and taking into account that the dominant behavior of K 3 will be K 3 ~ ksTCpPo, we 
remain with 

72 4 



F{d) ~ Ld 



d ' ^ id > (20) 



72^2 



\ k B T d(d 2 

The solution of this differential equation leads to the scaling form for the fluctuation free 
energy i.e. T(d) ~ d^ 2 ^ 3 . 



In the opposite limit we need to estimate K\ which is following de Gennes and Prost fl7 
obtained as K\ ~ ^ 3 ^-, where a is the hard core diameter of the polymers. Thus in this case 
we obtain 



T{d) ~ Ld 2 



(21) 



72^2 



\ k B T d(d 

the solution of which leads to the scaling form J-(d) ~ d~ 2 . By analyzing whether q± max 



is smaller or larger then 2.I^M^ it is easy to ascertain that the first limiting law should be 



valid for larger and the latter for smaller relative densities of the polymers. 

Interestingly we just derived the same scaling laws as those from corresponding mean 



field models by Helfrich [15] and De Gennes [16 . 



III. MEAN FIELD MODELS 

The mean field theory that corresponds to our problem describes a single polymer in a 
liquid crystalline matrix. In this view the polymer fluctuates in the field of its neighbors. 
These mean field models ( one could also say "Einstein" cage models ||21|| ) are strictly valid 
only for a hexagonal LC phase where the average position of a single polymer chain is 
well defined. In nematic-like phases the mean-square displacements are infinite and one 
has to choose a different language (eq. |l|). Although mean field models are not the proper 
description for our problem, they still illustrate the underlying physics. We therefore review 
their results. 

Let us imagine a flexible (elastic) polymer in a hard tube of diameter D, oriented on 



average along the axis of the tube, described with a Hamiltonian of a persistent chain [19 



H=\k B TC p $ds(^P ) \ (22) 

where C p is the persistence length of the polymer, associated with its bending elas- 
tic constant k c = k B TC p . For small deviations away from the average polymer posi- 
tion along the z-direction the two-dimensional displacement vector r(s) can be written as 
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r(s) = (r x (s),r y (s)). We assume that the diameter of the tube is small enough (smaller 
than the persistence length C p ) so that between two consecutive hits with the walls of the 
tube, the polymer propagates ballistically 

<{v{s)-v{s')) 2 >^C p {s-s'f. (23) 

The longitudinal correlation length Cn can be obtained from evaluating the angle of the hit 
with the wall 6 through statistics of a persistent chain i.e. 

D 2 £n 

9 2 ~ -2 ~ -1 wherefrom C\\ ~ £j /3 L> 2/3 . (24) 
|| p 

The free energy corresponding to the bumping between the polymer and the wall on a 
lengthscale defined by C\\ can thus be obtained as 

^^^£-2/3. (25) 



This result has been obtained previously by Odijk |20| and Helfrich Jl5| and represents a 
proper dimensional generalization of the Helfrich fluctuational force to 1-D confined elastic 
objects. 

On much larger length scales where the polymer chain effectively behaves as a free flight 
chain in 2-d, the effect of local nematic order can be modeled through a simple nematic 
coupling term in the elastic energy of the type (Nn(s)) 2 where N is the average nematic 
director and n(s) is the local polymer director, i.e. n(s) = r(s). If N is directed again along 
the z axis, the dominant part of the chain Hamiltonian then assumes the form 

K = h I *> (26) 

where g gives the strength of the nematic coupling. This Hamiltonian gives rise to a diffusive 
propagation of the polymer chain described with 

<(r( S )-r( S ')) 2 >^^ 1 (s-3'), (27) 

if compared to the ballistic propagation of a persistent chain Eq. |23|. For too large deviations 
the effective tube again simulates the effects of steric interactions. The angle of the hit with 
the wall of the tube is obtained as 
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D 2 1 

9 2 ~ — T ~ — — , wherefrom C\\ ~ p -D 2 , (28) 
£jj p£|l' II y ' v ; 

where we have now used the statistics of a free flight chain. The corresponding free energy 
is now obtained as 

JF ~ ^ ~ D -\ (29) 



Both results have already been derived for confined polymers ||16|| . The point we wanted 
to make is that they are consistent with our macroscopic theory if the cutoff is chosen 
appropriately, meaning q± max ~ c/ _1 while also equating the tube diameter D with the mean 
separation between the polymers d. 

A similar approach can be used also in the case where confinement is not described by 
a short-range hard-core interaction but is mediated via a soft confinement potential of the 
form V(r(z)). In this case, assuming that the fluctuations of the polymer position from its 
average are small, the Hamiltonian describing polymer elastic fluctuations can be written in 
the form 

1, r L , (d 2 v(z)\ 2 r L 



1, f L , fd 2 r(z)\ 2 1 „ r L 



- 2 H dz [^r) +1 2 v "i dzAz) (30) 

with 

^vi^)>=(i™ + ^)u, ( 31 ) 

where we have expanded the potential to second order in deviations from a straight line. 
One can define a length I* analogous to the Odijk length by minimizing the energy with 
respect to the longitudinal size £ of typical fluctuations |21j 



T ~ 7 (¥■ + v " rH ) — ci ~ (^) • (32) 



Thus we obtain for the confining energy 



L (V"\' 
F{D) ~ A;T— ~ k B TLy—J (33) 
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which apparently depends on the fourth root of the confining potential stiffness V" . This 



result is completely consistent with the one derived from a macroscopic theory eq. 15, 16 



in 



the limit q± ma x *C 2^^£, if we allow for the proper correspondence between the confining 
potential stiffness V" and the isothermal compressibility modulus B. 

If on the other hand the polymer is presumed to be completely flexible with an orien- 
tational part of the potential energy described with a nematic coupling of the form eq. |26|, 
the Hamiltonian of a single confined chain will contain a softer elastic part proportional to 
{~ds*~) ' I 11 ^ s case we have 

where the confining potential stiffness has been defined in the same manner as before. Just 
as before we can again introduce the appropriate Odijk length I* by minimizing the energy 
with respect to the longitudinal size I of typical fluctuations yielding p2| 

T ~iicr + v " rH ) — r2 ~(f)' < 35 > 

From here on by the same argument as before the confining energy scales as 



L [V"\ 1/2 
F(D)~kT-~k B TL( — \ (36) 



therefore as the square root of the confining potential stiffness. This result is now equivalent 
to the one derived from the macroscopic theory eq. [15], [16] but this time in the limit q± m ax ^ 



z y «i ■ 

To summarize, the mean field theories surprisingly reproduce the same scaling laws under 
similar conditions (high and low density regimes) as our more detailed calculation from part 
II. This indicates that the fluctuation part of the free energy is predominantely determined 
by the entropy loss of confining the polymer chains. Still, in order to model the free energy 
in more detail it is more appropriate to use our macroscopic model, because it predicts the 
correct crossovers between the different regimes. 

The mean-field approach with nematic coupling constant g as a free parameter, inferred 



from experiments, was used for the first time in |22[] to argue that the effect of elastic 
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fluctuations was to renormalize the decay length of underlying exponential repulsion. With 
the right choice of g, the calculated magnitude of the charge on DNA agrees reasonably 
with the numbers derived below and in ref. J23j. Odijk |24| also proposed what amounts to 



a variation on the mean-field theme. It is based on a variational estimate for the Gaussian 
width of a single chain density distribution function. This theory in general does not lead 
to a straightforward renormalization of the decay length of the underlying soft exponential 
interactions between chains. Also the phase boundaries calculated from this theory coupled 
to Lindemann's criterion |25] fall off the mark |26 ]. 



IV. MICROSCOPIC INTERACTIONS BETWEEN DNA MOLECULES 

In order to use our main result (eq. pl| |16|) for the equation of state of nematic polymer 
liquid crystals we first have to guess how the macroscopic elastic constants Ki,K 2 ,K 3 and 
lateral compressibility B depend on the density (po = 1/A = y/3/d 2 for hexagonally packed 
chains). Our first guess is the "one constant approximation" K\ = K 2 = K s « U(d)/d ||17|| , 
where U(d) is the interaction energy and d is the average distance between the molecules. 
In our case this approximation has to be modified, because the intrinsic bending stiffness k c 
of the DNA molecules contributes additionally to the bending Franks constant K 3 . 

K X = K 2 ^ U{d) jd 
K 3 ^ p k c + U(d)/d 

dV 2 " 4L \d 2 d d dd J y } 

We now ask about the intermolecular interactions between two DNA molecules. There 
are two major contributions to the repulsion between DNA molecules in monovalent salt 
solutions |p3]| : (1) screened electrostatic repulsion from negative charges along the DNA 
backbone; (2) hydration repulsion coming from partially ordered water close to the DNA 
surface p7|. At the ionic strengths and polyelectrolyte densities considered in this work 
there appears to be no important contribution to the attractive part of the total DNA-DNA 
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interaction: van der Waals forces are negligible |27] and counterion-correlation forces p8|j29 
are screened [[JIJ. 

The mean-field electrostatic interaction is best described by the Poisson - Boltzmann 
theory resulting in an interaction potential between two parallel charged rods ]31[] of the 
form 

U(d)/L = J?—K Q (d/\ D ) (38) 
2n€€o 

where is the Debye screening length. For ionic strengths / from monovalent salts, Xd = 



3.08A/J I[M]. Eq. |3B] refers to two infinitely thin line charges with a charge density per 
length £. This line charge density is related to the actual surface charge density a on a 
cylinder with radius a as follows 

£= 2na\ D /K 1 (a/\ D ) (39) 
For large d/X^ the Bessel function K can be approximated by 

9^777= W 

z - 'd/AD 



The hydration repulsion between solvated molecules in water can be described by the 
same formalism as for screened electrostatic repulsion (see The interaction energy 

therefore goes as K (d/\H) with X H m 3A p3f . 

Knowing about the intermolecular interaction we can now revisit our main result 



eq. p~5|J16| . From eq. [37] it is clear that Ki, K 2 and the bulk compressibility modulus perpen- 
dicular to the chains B essentially decay exponentially with d whereas K 3 after an initial 
exponential decay will at low densities be dominated by the intrinsic bending stiffness k c of 
the polymers. This fact has instructive consequences. For lower densities, where we would 



expect fluctuations to be more prominent, only the limiting form valid for q_L max ^ 2 /s^a 



remains. In this limit, using exponentially decaying interactions for B, the fluctuation part of 
the free energy goes essentially as the fourth root of the direct interaction, since K3 = pok c . 

Summarizing, the bare interaction can be described in the following way (for simplicity 
we use the free energy per length Q = T jV). 
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/=- -d/X H pfr p -d/Xn 

-±^= + bJ-^ 1 = (41) 

where a and b are the amplitudes of the hydration and screened electrostatic repulsion. The 
total free energy Q, using the low density limit eq. [L^ using q± max = n/d, is then 



riA\ r ia\ i u nr ( 7r ) 5/2 i -i/4 j ^go l^o ,„ x 
0(d) = ^o(d) + ^^-^ ^ w - (42) 



V. MATERIALS AND METHODS 

Sample preparation and determination of interaxial separations between DNA molecules 
by either x-ray or direct density measurements was described in detail previously ||. 

At very low osmotic pressures (1/100 atm) the elasticity of the dialysis bags could 
contribute to osmotic pressure that acts on the sample. For this reason we performed 
experiments in which we dialysed low concentration (0.1wt% — lwt%) Dextran solutions 
against each other, to make sure that there were no residual osmotic pressures resulting 
from partially inflated dialysis tubes. After equilibration the dextran concentrations inside 
and outside the dialysis bags agreed within 1% of the bathing concentration down to 0.1wt% 
Dextran. 

To compare our experimental results with our theory we expressed all data in terms of the 
interaxial spacing between two DNA molecules. Assuming hexagonal packing in all density 
regimes the relation between density and the interaxial spacing d is p = (610/d[A)) 2 ) [mg/ml]. 



The relation between the osmotic pressure II and Q is then ||27|| : 

% = ( 43 ) 
od 

The nonlinear fits were done using the Levenberg-Marquardt method, implemented in 
the data analysis software Igor 3.03 (WaveMetrics, OR). The fit function used was (using 
Qq from eq. f4"l"|) : 



eg og /4 d ja*g 1 ag (AA , 

dd = ^d + ckBTK V W - d W (44) 
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where a and b are the bare amplitudes of the hydration and the screened electrostatic 
repulsion. Since the prefactor of the fluctuation part of the free energy depends on the 
cutoff wave- vector we chose to fit it by the dimensionless constant c. The actual fit was 
performed using logio(dQ/dd) versus interaxial spacing d. We did that to achieve an equal 
weight of all data points over the whole osmotic pressure regime ( the osmotic pressure goes 
roughly exponentially in PEG-concentration ). 



VI. EXPERIMENTAL RESULTS 



Fig. 2 shows the equation of state (II — d) for long DNA molecules at different NaCl 
concentrations (lOmM - 2M). The figure presents a compendium of data (from [27j, [p3fl , f34| 



and heretofore unpublished data) obtained up to date on the upper portion of the DNA phase 
diagram. At high osmotic pressures (II) all ionic strengths merge into the same curve: an 
exponential decay with a decay length of about 3A. We attribute this behavior to structural 
forces in water (hydration forces) commonly observed between hydrated surfaces in water 



32]] . At lower osmotic pressures the curves start to deviate from each other reflecting 
the influence of screened electrostatic repulsion. Interestingly for ionic strengths > 1M 
the curves are independent of ionic strength over the whole osmotic pressure-regime. This 
indicates that for / > 1M the electrostatic contributions are sufficiently screened so that 
the equation of state is dominated by hydration repulsion alone. 

Fig. 3 shows the measured dQ/dd = \f?>Yid at 1M NaCl. Since from 1M on the 
electrostatic contribution is negligible we can use the data to determine the decay length of 
the hydration repulsion. We fitted the data according to eq. [44] using only one exponentially 
decaying direct interaction. As fit parameters we used the amplitude of the hydration 
repulsion a, the hydration decay length Xh and the prefactor c of the fluctuation part of the 
free energy. The fit is shown as a solid line and describes the data very well. The resulting 
hydration decay length was Xh = (2.9 ± 0.2) A. The amplitudes and the prefactor c are 
summarized in table 1. 
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Fig. 4 and Fig. 5 show measured dQ/dd for 0.5M and 0.1M NaCl. The solid line shown 



is the fit to eq. using a,b and c as fit parameter. The hydration decay length was set to 



\h = 2.88A and the Debye screening length was set to Ad = 3.08A/ y/[M]. The results are 
shown in table 1. 

Fig. 4 indicates the phase boundaries between the various liquid crystalline phase of 
DNA in detail. For all other salt concentrations we only indicated the isotropic to anisotropic 
transition. At very high osmotic pressures (not shown in Fig.4) there exists a crystalline 
hexagonal phase of DNA (see [37|). It melts into a line hexatic phase (regime a): a three 



dimensional liquid with long-range bond-orientational order perpendicular to the axis of 
the molecules |p8|| . As far as we can see, there is no indication for a hexagonal liquid 
crystalline phase in between the crystalline and the line hexatic phase. Between the line 
hexatic and the chiral phases of DNA the x-ray structure factor shows two peaks (regime 
b): a sharper peak (continuing regime a) at smaller interaxial spacings and a more diffuse 
one at wider spacings (continued by regime c). Since by using the osmotic stress method we 
hold all intensive variables (p,T,/j,) fixed the mesaured S(q) should originate from a single 
phase (Gibb's phase rule). At this point it is not clear whether the structure in regime b 
corresponds to a new phase in between a line hexatic phase (nonchiral) and a cholesteric 
phase. 

The chiral phases start with a cholesteric phase as observed by electron microscopy 
and polarization microscopy fl39l,f40f . At even lower concentrations there is some indication 



for additional chiral phases (precholesteric phase |^T| and blue phases [f|2]]). Finally the 



anisotropic liquid crystalline phase melts into an isotropic phase |26j (regime d). 

VII. DISCUSSION 

DNA is highly charged ( two negative charges per base-pair or 3AA of its length). In an 
electrolyte solution the DNAs net negative charge creates an accumulation of counterions 
close to its surface that screen part of the bare charge and lead to an effective charge density 



18 



that is felt at long distances between chains. Theoretically this effect can be captured by 
nonlinear Poisson-Boltzmann theory |53| . 



Even though our nematic polymer liquid crystalline theory was based on a rather simplis- 
tic model it describes all the data fairly well with reasonable values for the fitted parameters. 
In table 1 we have summarized the corresponding line-charge densities and surface-charge- 



densities that were calculated according to eq. |38] and eq. |39|. At 0.5M and 0.1M the fits 
give the same surface-charge density a = 0.07 C/m 2 . This value corresponds to about 50% 
of the bare charge of DNA (0.15 C/m 2 ). 

Our result of 50% effective charge agrees very well with the analysis of electrophoretic 
measurements jPj by Schellman and Stigter |45| that resulted in about 60% effective charge 



density for Na-DNA. A recent study using steady-state-electrophoresis reported 10% for Na- 
DNA ||46|| . This discrepancy exists not because of experimental uncertaintly, but because of 
different theoretical treatments of the measured mobilities (for a recent discussion see [§j]). 
Part of the difference might come from the complicated details of frictional forces acting on 
a rough, charged polyelectrolyte. 

One has to bear in mind that for the 0.5M data the uncertainty of the fitted amplitude 
is almost 100%. This uncertainty is not very surprising considering the close proximity of 
the two decay lengths Xh = 2.9 A and Xd = 4.36A. As soon as the decay lengths separate 
from each other, as in the case for 0.1M, the statistical uncertainty drops to 10%. 

We also fitted the prefactor c. If in the theory (eq. ^) one chooses the cutoff wave- vector 
to be at the Brillouin-zone radius g_i_ max = tr/d, c evaluates to 0.16. The fitted pref actors 
c range from 1.3 to 0.8. This is about 5-8 times larger than the theory predicted. On the 
other hand, considering the simple underlying model, the value is not too far off. Choosing a 
cutoff at twice the Brillouin-zone radius, for example, would give a value right on the fitted 
one. In our view, the prefactor depends on the fine details of the in-plane structure-factor 
S(q±) (see eq. [19]). In the case of exponential direct interactions between the chains any 
algebraic dependence of the cutoff with respect to the density will be dominated by the 



fourth root of the direct interaction (eq. |42|). The fact that for different ionic strength the 
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pref actors have similar values ( all around 1 ) strengthens our argument. 

The fits even seem to prove explicitly the presence of hydration repulsion. If one tries to 
calculate the charge densities for 1M, using a screening length of Ad = 3.08A and ignoring 
any contribution of hydration repulsion, the surface charge density results in a = 0.19 C/m 2 , 
25% more charge than the total phosphorus charges on DNA. Since it is generally believed 
that more than half of the bare charge on DNA is screened, a picture with pure electro- 



static double layer repulsion is hard to envision. Recently Lyubartsev and Nordenskiold f48 
published Monte-Carlo simulation addressing the DNA case, comparing their results with 
osmotic stress measurements from our lab. They showed that electrostatic repulsion strongly 
increased beyond simple linearized PB theory, when two charged cylinders approach each 
other closely. However this simulation used a solid cylinder-model of DNA and hard sphere 
d = 4A ions. The measurements are in a regime where most of the aqueous volume is inside 
the DNA grooves rather than outside any DNA cylinder. Better models have to include the 
possibility of ions entering the groove space. Although their results seem to fit the data 
for 0.5M quite well, there is not as good success at other salt concentrations. We take the 
force data in high salt concentration to be strong evidence for hydration repulsion. Even 
non-charged polymers, like the polysaccharide schizophyllan show an exponentially decaying 
repulsion with a decay length of about Xh = 3.4A [fj9"| . Hydration repulsion is a general 
feature of water soluble molecules at separations < lnm |32] that can not be simply denied 

sg. 



Osmotic stress measurements can be used to determine more directly effective charge 
densities of semiflexible polyelectrolytes. We can therefore test and compare theories that 
predict effective charge densities, like non-linear Poisson-Boltzmann (PB) or Manning theory 
f5T| . Our results indicate that for Na the effective charge is twice as large as predicted by both 



theories. Previous studies in our lab |23|] observed significant differences in charge densities 
using different counterions, such as Li,Na,K,Cs,Tri-methyl-ammonium. These results merit 
further analysis. 

The form of the fluctuation part of the free energy (eq. ^2|) suggest that fluctuation 
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enhanced repulsion may be important for many lyotropic polymer liquid crystalline systems. 
The prefactor c of the fluctuation part only depends weakly on the bending constant of the 
polymer c oc k~ l / A . The only condition for enhanced repulsion is that the polymers are 
longer than their persistence length, and that they remain in a nematic-like phase. 

Another appealing conclusion from our work is the emerging connection between multi- 
lamellar lipid and columnar polyelectrolyte (DNA) arrays. Both are governed by the same 
type of colloidal forces |52 ] , except that for polyelectrolyte arrays the attractive forces at rel- 



evant ionic strengths and polymer concentrations are usually negligible, and show the same 
conformational flexibility describable by an elastic term in the conformational Hamiltonian 
(Helfrich ]53| in the case of lipid multilayers and Kratky-Porod [19 in the case of semiflexible 



chains). The interplay between fluctuations and effective interaction in multilamellar arrays 



has been a topic for quite awhile |54| and is at present reasonably well understood. 

Setting aside the attractive part in the bare interaction potential and the different di- 
mensionality of the fundamental interacting objects (2-D in the case of lipid multilayers 
vs. 1-D in the case of polyelectrolytes) , multilamellar lipid and columnar polyelectrolyte 
arrays share exactly the same Hamiltonian. This leads to the same type of fluctuational 
renormalization of interactions force that can be expressed by the following two forms of 
free energy 

k B TV r& 5/2 , _ k B TV [E , 




•^5x23/% fe 9lmai + '" VS ' J "l67te 9lmm + -' (45) 



The first one pertains to the fluctuations in polyelectrolyte arrays, c.f. eq. [16|, and the second 
one to the fluctuations in multilamellar systems, e.g. ||55|| . Both these results are written in 
the limit of small density of fluctuating objects. The apparent formal differences between the 
two expressions are due solely to the dimensional difference between the fluctuating objects, 
i.e. 1-D as opposed to 2-D. 

The main difference between multilamellar and columnar arrays interacting through 
exponential repulsive forces would thus be the two-fold vs. four-fold renormalisation of the 
decay length. 
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There have recently (see [56[ for an overview), and maybe not so recently j57|, appeared 
quite a few speculations on the possible attractive component to the polyelectrolyte inter- 
action forces. At polyelectrolye densities and ionic strengths described in this work, there 
is certainly no evidence to presume there are any. We can not, however, exclude the possi- 
bility that low ionic strengths and low polyelectrolyte densities, thus promoting pronounced 
unscreened counterion fluctuations [p8j^0| , conspire to bring forth non-negligible attractive 
forces between DNA molecules. Should this turn out to be the case, columnar polyelectrolyte 
arrays would become even more similar to multilamellar lipid arrays. DNA and cell mem- 
branes are among the principal organizational structures in biology. That they share such 
a pronounced amount of common physics certainly add up to a rather pleasing intellectual 
development. 
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FIGURES 

FIG. 1. Illustration of how splay couples to the density of different nematic liquid crystals: a) 
stiff, long rods; b) stiff, short rods; c) semiflexible, long polymers. Long, stiff rods (a) show strong 
density changes when splayed. In the case of short and stiff rods (b) the voids created by splaying 
the material are filled by other short rods. For long, semiflexible polymers (c) the voids are filled 
by polymers folding back on itself. 

FIG. 2. Equation of state for DNA liquid crystals at ionic strength from 150mM to 2M NaCl. 
We plotted logli versus interaxial spacing d. The interaxial spacing d was measured by x-ray 
scattering. 

FIG. 3. Measured dQ/dd for DNA liquid crystals at ionic strength of 1M over 5 orders of 
magnitude in osmotic pressure. The solid curve represents a fit of the datapoints in the anisotropic 
regime to a nematic liquid crystalline theory (see text) assuming exponential repulsion (hydration 
and screened electrostatics). The broken line represents the bare interaction without fluctuation 
enhanced repulsion. Interaxial spacings above the -.- line are measured by x-ray scattering; below 
this line the spacing d are derived from measured DNA densities where molecules are expected to 
be hexagonally packed. 

FIG. 4. Measured dQ/dd for DNA liquid crystals at ionic strength of 0.5M over 5 orders of 
magnitude in osmotic pressure. Four structural regimes (a-d) could be distinguished. See Fig.3 for 
annotations. 

FIG. 5. Measured dQ/dd for DNA liquid crystals at ionic strength of 0.1M. See Fig.3 for 
annotations. 
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TABLES 

TABLE I. summary of all fitted parameters for 1M, 0.5M and 0.1M NaCl, as well as the corre- 
sponding line charge density £, surface charge density a and the fraction of unscreened phosphorus 
charges E 



I[M] 


a[J/m] 


Xh [A] 


b[J/m] 


Ad [A] 


c 


S[C/m\ 


a[C/m 2 } 


E 


1 


(1.7±0.9)10~ 7 


(2.9 ±0.2) 




3.08 


(1.2 ±0.1) 








0.5 


(1.4±0.5)10- 7 


2.88 


(3±2.6)10- 9 


4.36 


(1.3 ±0.2) 


2.1 • 10~ 9 


0.07 


0.49 


0.1 


(1.1±0.3)10~ 7 


2.88 


(4.1±0.3)10~ 10 


9.74 


(0.8 ±0.06) 


7.8- 10~ 10 


0.07 


0.48 
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